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Abstract 

Di-jet production is studied in collisions of quasi-real photons radiated by the LEP beams at 
e’''e“ centre-of-mass energies = 161 and 172 GeV. The jets are reconstructed using a cone 
jet finding algorithm. The angular distributions of direct and double-resolved processes are 
measured and compared to the predictions of leading order and next-to-leading order pertur¬ 
bative QCD. The jet energy prohles are also studied. The inclusive two-jet cross-section is 
measured as a function of and and compared to next-to-leading order perturbative 
QCD calculations. The inclusive two-jet cross-section as a function of is compared to the 
prediction of the leading order Monte Carlo generators PYTHIA and PHOJET. The Monte 
Carlo predictions are calculated with different parametrisations of the parton distributions of 
the photon. The influence of the ‘underlying event’ has been studied to reduce the model 
dependence of the predicted jet cross-sections from the Monte Carlo generators. 
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1 Introduction 


We present a study of di-jet production in photon-photon collisions at -y/s^g = 161 and 172 GeV 
with an integrated luminosity of 20 pb“^. The cone jet finding algorithm was used to reconstruct 
jets. The production of di-jet events in the collision of two quasi-real photons can be used to 
study the structure of the photon and to test QGD predictions. At e'''e“ colliders the photons 
are emitted by the beam electron^. Most of the photons carry only a small squared four- 
momentum, and can be considered to be quasi-real ~ 0). Accordingly, the electrons 

are scattered with very small angles and are not detected. Events where one or both scattered 
electrons are detected are vetoed (“anti-tagged”). 

The interactions of the photons can be modelled by assuming that each photon can either 
interact directly or appear resolved through its fluctuations into hadronic components. The 
interaction of two photons can be classihed either as a direct process where two bare photons 
interact, a single-resolved process where a bare photon interacts with a parton (quark or gluon) 
of the other photon or a double-resolved process where partons of both photons interact. The 
possibility to distinguish between direct and resolved processes in di-jet events has already been 
demonstrated by OPAL at y^gg = 130 — 136 GeV |^. Depending on the type of photon-photon 
interaction different matrix elements for the QGD scattering process contribute. These matrix 
elements have been calculated in leading order (LO) |0] and next-to-leading order (NLO) 0. 
The matrix element of the scattering between two bare photons is the one for the process 
77 9.Q- III double-resolved processes the matrix elements of quark-quark, gluon-quark and 

gluon-gluon scattering are involved 0 . These calculations predict different distributions of the 
parton scattering angle in the centre-of-mass system of the colliding particles. In hadron-hadron 
and photon-hadron interactions similar QGD predictions have already been conhrmed |^,Q. 

The investigation of the internal structure of jets gives insight into the transition between 
a parton produced in a hard process and the observable hadrons which originate from the 
fragmentation process . The dependence of the jet shapes on QGD parton radiation calculated 
in the Leading-Log Approximation (LLA), and the differences between the jet shapes of direct 
and resolved processes, have recently been measured at HERA [0. 

1 Positrons are also referred to as electrons 
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The measurement of inclusive jet cross-sections in 77 and 7 p interactions can constrain 
the gluonic content of the photon This is done by comparing the jet cross-sections to 

the LO QCD models PYTHIA and PHOJET, using different parametrisations of the parton 
distribution functions of the photon, and to NLO QCD calculations |I0|, |Tl|] . PYTHIA and 
PHOJET also model the so-called ‘underlying event’ by multiple interactions involving partons 
from the remnants of the same two initial photons, whereas the NLO QCD calculations do not 
take into account such effects. In the models, the contribution from multiple interactions to 
the jet cross-sections has to be tuned to increase sensitivity to the parton distributions of the 
photon. In contrast to deep inelastic electron-photon scattering |^, which in leading order 
is only sensitive to the quark content of the photon, the gluon content of the photon can be 
tested directly in the resolved interaction of two almost real photons, where the fraction of 
gluon initiated processes is large. 


Inclusive one-jet and two-jet cross-sections in photon-photon collisions have previously been 
measured at an e’''e“ centre-of-mass energy of -y/Sg^ = 58 GeV at TRISTAN ||,|^ and at an 
e+e“ centre-of-mass energy of = 130 and 136 GeV at LEP |l|. This paper extends our 
analysis at lower energies in which a similar strategy was used. Jets are studied in a wider 
kinematic range and with higher integrated luminosity. In addition, we present new results on 
angular distributions, jet shapes and energy flows. 


2 Process kinematics 


The properties of the interacting photons are described by their negative squared four-mo¬ 
mentum transfers, Qj. Each is related to the electron scattering angle 0' relative to the 
beam direction by 

~ - cose'), (1) 

where pi and p' are the four-momenta of the beam electrons and the scattered electrons, respec¬ 
tively, and Ei and E' are their energies. Events with detected scattered electrons (single-tagged 
or double-tagged events) are excluded from the analysis. This anti-tagging condition dehnes 
an upper limit on for both photons. This condition is met when the scattering angle 0' of 
the electrons is less than 0 max) where 0 max angle between the beam-axis and the inner 

edge of the detector. The squared invariant mass of the hadronic hnal state 

-fep-) (2) 


is calculated by summing over the energies. Eh, and momenta, ph, of all hnal state hadrons. 
The spectrum of photons with an energy fraction y of the electron beam may be obtained by 
the Equivalent Photon Approximation (EPA) |I^: 


f'y/e{y) 


^ ( 1 + {l-y? , ^ 

271 [ y ^ QL 
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with a being the electromagnetic coupling constant. The minimum kinematically allowed 
squared four-momentum transfer, is determined by the electron mass me- 

2 ^ 

Vmin 1 

1-2/ 
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The effective maximum four-momentum transfer given by the anti-tagging condition, 

i.e. the requirement that both electrons remain undetected. 


3 The OPAL detector 


A detailed description of the OPAL detector can be found in Ref. [15], and therefore only a 


brief account of the main features relevant to the present analysis will be given here. 


The central tracking system is located inside a solenoidal magnet which provides a uniform 
axial magnetic held of 0.435 T along the beam axisQ. The magnet is surrounded in the barrel 
region (| cos6\ < 0.82) by a lead glass electromagnetic calorimeter (ECAL) and a hadronic sam¬ 
pling calorimeter (HCAL). Outside the HCAL, the detector is surrounded by muon chambers. 
There are similar layers of detectors in the endcaps (0.82 < | cos6\ < 0.98). The small angle re¬ 
gion from 47 to 140 mrad around the beam pipe on both sides of the interaction point is covered 
by the forward calorimeters (FD) and the region from 25 to 59 mrad by the silicon tungsten 
luminometers (SW) JT^. From 1996 onwards, relevant to the data presented in this paper, the 
lower boundary of the acceptance has been increased to 33 mrad following the installation of a 
low angle shield to protect the central tracking system against possible synchrotron radiation. 


Starting with the innermost components, the tracking system consists of a high precision 
silicon microvertex detector, a vertex drift chamber, a large volume jet chamber with 159 
layers of axial anode wires and a set of ^ chambers measuring the track coordinates along 
the beam direction. The transverse momenta, p^, of tracks are measured with a precision 
parametrised by a^^jp^ = y'^ 0.022 -|- (0.0015 ■ px)^ (pt in GeV/c) in the central region. In this 
paper “transverse” is always defined with respect to the ^ axis. The jet chamber also provides 
measurements of the energy loss, dE/dx, which are used for particle identification [1^ . 


The barrel and endcap sections of the ECAL are both constructed from lead glass blocks 
with a depth of 24.6 radiation lengths in the barrel region and more than 22 radiation lengths 
in the endcaps. The FD consist of cylindrical lead-scintillator calorimeters with a depth of 24 
radiation lengths divided azimuthally into 16 segments. The electromagnetic energy resolution 
is about 18%/vT, where E is in GeV. The SW detectors consist of 19 layers of silicon detectors 
and 18 layers of tungsten, corresponding to a total of 22 radiation lengths. Each silicon layer 
consists of 16 wedge shaped silicon detectors. The electromagnetic energy resolution is about 
25%//F {E in GeV). 


4 Event selection and jet finding 


Two-photon events are selected with the following set of cuts: 


• The sum of all energy deposits in the EGAL and the HGAL has to be less than 45 GeV. 

^In the OPAL coordinate system the z axis points in the direction of the e“ beam. The polar angle 0, the 
azimuthal angle (j) and the radius r denote the usual spherical coordinates. 
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Calorimeter clusters have to pass an energy threshold of 100 MeV in the barrel section or 
250 MeV in the endcap section for the ECAL and of 600 MeV for the barrel and endcap 
section of the HCAL. 

• The visible invariant mass measured in the ECAL has to be greater than 3 GeV. 

• The missing transverse energy of the event measured in the ECAL and the ED has to be 
less than 5 GeV. For a ED cluster to be counted its energy has to be larger than 1 GeV. 

• At least 5 tracks must have been found in the tracking chambers. A track is required 
to have a minimum transverse momentum of 120 MeV/c, at least 20 hits in the central 
jet chamber, and the innermost hit of the track must be within a radius of 60 cm with 
respect to the z axis. The distance of the point of closest approach to the origin in the 
rcj) plane must be less than 30 cm in the z direction and less than 2 cm in the ref) plane. 
Tracks with a momentum error larger than the momentum itself are rejected if they have 
fewer than 80 hits. The number of measured hits in the jet chamber must be more than 
half of the number of possible hits, where the number of possible hits is calculated from 
the polar angle 9 of the track, assuming that the track has no curvature. 

• To remove events with scattered electrons in the ED or in the SW calorimeters, the 
total energy sum measured in the ED has to be less than 50 GeV and the total energy 
sum measured in the SW calorimeter has to be less than 35 GeV. A cluster in the SW 
calorimeter is accepted if it has an energy of more than 1 GeV. These cuts also reduce 
the contamination from multihadronic annihilation events with their thrust axis close to 
the beam direction. 

• In order to estimate the position of the primary vertex, we calculate the error-weighted 
average {zq) of the ^ coordinates of all tracks at the point of closest approach to the origin 
in the r(f> plane. The background due to beam-gas interactions is reduced by requiring 
I ( 2 : 0 ) I < 10 cm and \Q\ < 3, where Q is the net charge of an event calculated from adding 
the charges of all tracks. 

• To remove beam-wall events the radial distance of the primary vertex from the beam axis 
has to be less than 3 cm. 


In the cone jet finding algorithm the total transverse energy of the jet inside the cone 
is the scalar sum of the transverse energies of its components |l|, |^. In all parts of this 


analysis, a sum over the particles in the event or in a jet means a sum over tracks satisfying the 
above quality cuts, and over all calorimeters clusters, including the ED and SW calorimeters. 
An algorithm is applied to avoid double-counting of particle momenta in the central tracking 
system and their energy deposits in the calorimeters |l|. The transverse energy Et. of a 
particle i is dehned relative to the 2 ; axis of the detector with = Ejsin^j. For a cone jet 
to be accepted, the value of must be greater than a certain minimum E™™. The results 


of the cone jet Ending algorithm depend on and the cone size R = (Ar^)^ -|- (A^)^ with 

pseudorapidity rj = — lntan(6*/2) and azimuthal angle 0. Here the values were chosen to be 
R = 1 and = 2 GeV. The jet pseudorapidity in the laboratory frame is required to be 

within < 2. Monte Garlo studies have shown that jets with < 2 are well reconstructed 
and are normally fully contained in the detector. This leads to a wider acceptance than in 
our previous analysis |^. 
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We use data corresponding to an integrated luminosity of 9.9 pb ^ at = 161 GeV 
and 10.0 pb“^ at \/s^^ = 172 GeV. After applying all cuts and requiring at least two jets with 
> 3 GeV and < 2, 2845 events remain, equally divided between the two centre- 
of-mass energies. For the purpose of this analysis, the difference between the data taken at 
\/s^^ =161 GeV and at 172 GeV is small and therefore the distributions for both energies have 
been added. About 12.3% of the di-jet events contain 3 or more jets. In events with more than 
two jets, only the two jets with the highest 77^* values are taken. 


5 Monte Carlo simulation 


The Monte Garlo generators PYTHIA 5.722 |]^|^ and PHOJET 1.05c are used, both 
based on LO QGD calculations. These generators have been optimised to describe yp and pp 
interactions. The probability of Ending a parton in the photon is taken from parametrisations 
of the parton distribution functions. The SaS-lD parametrisation is used as default in 


PYTHIA and the LO GRV parametrisation as default in PHOJET. All possible hard 
interactions of quarks, gluons and photons are simulated using LO matrix elements for massless 
quarks. The default value of the cutofi on the transverse momentum of the two outgoing partons 
of p™™ = 1.4 GeV/c is used. More details can be found in Ref. |^. 


The incoming photons in double-resolved events can be viewed as beams of partons. For 
small parton transverse momenta the LO parton scattering cross-section diverges and becomes 
larger than the non-diEractive cross-section as measured in yp collisions. If more than one 
parton scattering process is allowed in one event the problem of too large parton cross-section 
can be solved. These multiple interactions are calculated as LO QGD processes between partons 
of the photon remnants. In PYTHIA a lower cutofi parameter is introduced, which describes 
the transverse momentum of the parton and is set by default to 1.4 GeV/c. The fluctuations of 
the number of hard interactions are calculated from a Poisson distribution. The average number 
Umi of interactions in double-resolved di-jet events simulated by PYTHIA is 1.3 for the default 
setting p™ = 1.4 GeV/c using SaS-lD. PYTHIA and PHOJET use multiple interactions as a 
component to model the underlying event. 


In PHOJET the suppression of the total yy cross-section is parametrised using Gener¬ 
alised Vector Meson Dominance (GVMD). A model for the change of soft hadron production 
and diffraction with increasing photon virtuality is also included. The photon-photon mode 
of PYTHIA only simulates the interactions of real photons with = 0. The virtuality of 
the photons defined by enters only through the Equivalent Photon Approximation in the 
generation of the photon energy spectrum, but the electrons are scattered at zero angle. This 
model is not expected to be correct for larger values of Q^. The contribution of di-jet events 
with > 1 GeV^ and 0' < 33 mrad generated with the electron-photon mode of PYTHIA is 
small and therefore neglected. 


The fragmentation of the parton final state is handled in both generators by the routines 
of JETSET 7.408 |]^. Initial and final state parton radiation is included based on the LLA. 


All signal and background Monte Garlo samples apart from beam-gas and beam-wall events 
were generated with full simulation of the OPAL detector [^|. They are analysed using the 
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same reconstruction algorithms as are applied to the data. 


The median of the selected PHOJET events is of the order 10“^ GeV^. The visible 
hadronic invariant mass, ITvis, measured with all detector information is well described within 
the errors of the measurements by the Monte Carlo simulations. A detailed comparison between 
Wvis and the generated W can be found in Ref. p3 . 


After applying the detector simulation and the selection cuts to these events, about 83% of 
all generated Monte Carlo events with at least two hadron jets in the range > 3 GeV and 
Ir/j®*! < 2 are selected. The trigger efficiency for all selected Monte Carlo events with at least 
two reconstructed jets in the detector is close to 100%. The number of background events is 
small, about 1.5% in total. About 1 . 1 % of the events in the data sample are expected to be 
e'''e“ annihilation events with hadronic hnal states and 0.4% e’''e“ —> e'^e~T~^T~ events. No 
signihcant background from beam-gas or beam-wall events is observed. 


6 Properties of direct and resolved processes 


In LO QCD, neglecting multiple parton interactions, two hard parton jets are produced in 77 
interactions. In single- or double-resolved interactions, the two hard parton jets are expected 
to be accompanied by one or two remnant jets. 


A pair of variables, and , can be dehned which specify the fraction of the photon’s 
momentum participating in the hard scattering: 


i: (E+p.) 

_l_ jets=l,2 , 

= V- /rn -5 ^7 

E (E + p,) 

hadrons 


i: (E-p,) 

jets=l,2 

i: (E-p.y 

hadrons 


( 3 ) 


where pz is the momentum component along the z axis of the detector and E is the energy of 
the jets or hadrons. These variables give some separation between direct and resolved di-jet 
events |II]. 


Ideally, for direct events without remnant jets, the total energy of the event is contained in 
the two jets, i.e. = 1 and x~ = 1, whereas for single-resolved events either x^ or x~ and 
for double-resolved events both values, and x~, are expected to be smaller than 1 . Samples 
with large direct and double-resolved contributions can be separated by requiring both and 
x~ to be larger than 0.8 (denoted as x^ > 0 . 8 ) or both values to be smaller than 0.8 (denoted 
as x^ < 0.8), respectively. Details about the separation between different event classes can be 
found in Ref. |^. In the PYTHIA Monte Carlo using the SaS-lD parametrisation 86 % of all 
events in the region x^ > 0.8 originate from direct interactions and 81% of all events in the 
region < 0.8 originate from double-resolved interactions. 


The Xy distribution is shown in Figure |l] in bins of where 
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is the mean value of the transverse energies and of the two jets. Each event is added 
to the plot twice, at the values of and of x~. No correction for selection cuts and detector 
effects has been applied, but the background has been subtracted using the Monte Carlo. The 
Monte Carlo predictions of PYTHIA and PHOJET are normalised to the number of events 
observed in the data. The contribution from direct processes, as predicted from PYTHIA, is 
also shown. The events from direct processes are concentrated at high values. In Fig. 0 a, 
at low the direct part contributes to about 17% to the total number of events. As Ep* 
increases, the x^ distribution shifts to higher values and the fraction of direct events in the 
PYTHIA sample increases to 68 % for 12 < < 20 GeV (Fig. |^). The number of events is 

underestimated by PYTHIA and PHOJET by about 25 — 30%, if the predicted Monte Carlo 
cross-sections are taken into account, mainly for x^ < 0.9. 


7 Angular distributions in direct and resolved events 


Since the jets in double-resolved events do not contain all of the hadronic activity it is ex¬ 
pected that there will be more energy flow outside the jets in double-resolved events than in 
direct events. Fig. § shows the energy flow transverse to the beam direction as a function of 
Ap' measured with respect to the jet direction for data samples with different x"^ cuts. The 
pseudorapidity difference is dehned by: 

Ap' = k{T] — 

where r] is the pseudorapidity of the cluster or the track. The factor k is chosen event-by- 
event to be A: = J-l for events with > x~ and fc = — 1 for events with x'^ < x~. As a 
consequence, there is always more of the remnant found at At]' < 0 and the enhancement due 
to the additional transverse energy flow observed at negative and positive Ap' is asymmetric. 
No correction for acceptance or resolution effects has been applied. The energy flow is integrated 
over |A0| < 7 r/ 2 . The jets in the data sample with > 0.8 (Fig. are more collimated and 
there is almost no activity outside the jet, whereas the transverse energy flow of two-jet events 
with < 0.8 (Fig. |]b) shows considerable activity outside the jets for jAr^'l > 1. The energy 
flow outside the jets is well modelled by PYTHIA whereas PHOJET shows a wider distribution. 

In the di-jet centre-of-mass frame one expects different angular distributions for direct and 
double-resolved events. An estimator of the angle 9* between the jet axis and the axis of the 
incoming partons or direct photons in the di-jet centre-of-mass frame can be formed from the 
jet pseudorapidities: 

( jetl _ jet2 

— 

2 

The two jets are assumed to be collinear in 0 with . Since the ordering of the jets in 

the detector is arbitrary, only | cos0*| can be measured. The matrix elements of elastic parton- 
parton scattering processes are known to LO p[. For a given parton centre-of-mass energy 
the cross-sections vary only with the scattering angle 9*. The LO direct process 77 — >■ qq is 
mediated by f-channel spin-| quark exchange which leads to an angular dependence oc (1 — 
cos^0*)“^. In double-resolved processes all matrix elements involving quarks and gluons have 
to be taken into account, with a large contribution from spin-0 gluon exchange. After adding 
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up all relevant processes, perturbative QCD predicts an angular dependence of approximately 
OC (1 — Icos^*!)”^ p. 


In order to measure the cos 6* distribution, additional cuts have to be applied. These cuts 
minimise kinematic biases and improve the detector resolution on cos 9*. The invariant mass 
of the two-jet system is calculated as 


Mjj 


2E-. 


ijet 


1 — Icos^* 


The cut on > 3 GeV restricts the accessible range of values of | cos6**|. For Mjj < 10 GeV 
the number of events decreases because of the cut. Requiring Mjj to be larger than 12 GeV 
ensures that the distribution in the range | cos6**| < 0.85 is not biased by the cut. The 
Lorentz boost of the two-jet system in the 2 ; direction is dehned by 


^jetl _j_ ^jet2 



since in the two-jet centre-of-mass system rj* = 0. Events with |rj| > 1 were rejected because 
the detector resolution on |cos0*| deteriorates signihcantly for events with \f]\ > 1. After 
additionally requiring Mjj > 12 GeV and \f]\ < 1 150 data events remain with > 0.8 and 
350 data events with < 0.8. 

Table |l| shows the measured cross-sections for data samples with large direct and large 
double-resolved contributions according to the separation with x^ and x~. The cross-section in 
each bin of | cos 6* \ was corrected according to the efficiency found for that bin in the PYTHIA 
and PHOJET samples. The central value is the mean of the result from PYTHIA and PHOJET. 
The systematic uncertainty on the jet cross-sections given in this section is determined by 
varying the energy scale of the ECAL in the Monte Garlo simulation by ±5%. The dependence 
on the Monte Garlo models used is taken into account by adding the difference between the 
results obtained with PYTHIA and PHOJET to the systematic error. The contributions of 
these two errors to the total systematic error are of similar size. The systematic errors due to 
the luminosity measurement, the trigger efficiency and the hnite number of Monte Garlo events 
are small in comparison. 


Fig. g shows the bin-by-bin corrected | cos 0*1 distribution of events with x^ > 0.8 and of 
events with < 0.8. The abscissae of the data points are plotted according to the method 


proposed in Ref. [26|. The predictions of the theoretical parton distributions are integrated to 


hnd the position of the data points. The error of this position is obtained using the predictions 
of the different parton processes. It is smaller than the line width. The predicted curves 
have been normalised to the data in the hrst three bins, in order to compare the shape of the 
measured cross-sections as a function of | cos 6* \ with the QGD matrix element calculation and 
the NLO QGD calculation The error bars show the statistical and the systematic errors added 
in quadrature. The overall error on the normalisation of the data points is dominated by the 
statistical error in the hrst three bins, which is about 20%. The events with x. 
a small rise with I cos 9* 


^ > 0.8 show 

whereas the events with x^ < 0.8 show a much stronger rise with 
cos 0*1, as expected from the QGD calculations. In the Monte Garlo events about 10% of the 


processes with x^ > 0.8 are double-resolved and about 15% of the processes with x^ < 0.8 are 
direct. The | cos0*| distribution is not much affected by these impurities, since the Monte Garlo 
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double-resolved events with > 0.8 show a smaller rise with | cos 6 **| than the double-resolved 
events with < 0.8 and the Monte Carlo direct events with x^ < 0.8 show a stronger rise 
with I cos 0*1 than the direct events with x^ > 0 . 8 . 

In Fig. ^ the points for < 0.8 lie close to the predictions of a QCD matrix element 
calculation of the interaction of quarks or gluons in the photon [Q . The matrix elements with a 
relevant contribution to the cross-section where anti-quarks are involved instead of quarks show 
a similar behaviour to the examples shown. The points for > 0.8 are comparable with the 
results of a calculation of the process 77 —qq. The QCD matrix element calculations agree 
well with the data samples with large direct and large double-resolved contribution. 

The data points are compared in Figure ^b to a NLO perturbative QCD calculation |Q. The 
contribution of the different processes to all double-resolved events depends on the parametrisa- 
tion of the parton distribution functions. This calculation uses the NLO GRV parametrisation 
and was repeated for the kinematic conditions of this analysis. The shape of the data points 
and the NLO calculation agrees well. However, before normalisation the predicted cross-section 
is a factor of two too high for the direct events and about 50% too small for double-resolved 
events. A NLO QCD calculation using the GS photon structure function |2^ (not shown) 
shows a similar behaviour. 


8 Jet shapes 


The internal structure of jets produced in photon-photon interactions is studied at the hadron 
level. The jet shape is characterised by the fraction of a jet’s transverse energy 
inside an inner cone of radius r concentric with the jet dehning cone [^]: 


'0(r) 


1 ^ ET^jr) 


(4) 


where E^{r) is the transverse energy within the inner cone of radius r and Ajet is the total 
number of jets in the sample. By dehnition, 'ip{r = R) = 1. 


To check the relative importance of parton radiation and fragmentation in the formation of 
a jet the parton-shower in the LLA as implemented in PYTHIA 5.7 |]^ has been used. The 
jet shape is affected both by fragmentation and gluon radiation. 


The jet shapes are corrected to the hadron level using the Monte Carlo event samples with 
single-resolved, double-resolved and direct processes. The corrected jet shapes are denoted by 
'ip{r) and refer to jets at the hadron level with a cone radius R = 1. The reconstructed jet 
shapes are corrected for acceptance effects and the hnite detector resolution. The correction 
factors also take into account the selection criteria and the purity and efficiency of the jet 
reconstruction. The corrected jet shapes are determined bin-by-bin as 


^lJ{r) = C'(r) ■ V’det(r), 


where the correction factors C{r) are dehned as 


C(r) = 
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Since correction factors are different for single-resolved, double-resolved and direct processes, 
the correction factors have been calculated for each process separately in each interval. For 
direct processes C{r) is close to 1 for all r, whereas for single- and double-resolved processes 
C{r) lies in the range from 0.8 to 1.2 for r = 0.2 and approaches 1 as r increases. The relative 
contributions of the different processes have been estimated by comparing the distribution 
of the data events with PYTHIA (see Figure |l|) and changing the contributions to get a correct 
description. In the lowest bin (3 < < 6 GeV) PYTHIA describes the x.^ distribution 

well and the contributions of the different processes were not modihed. For the bins in the 
range 6 < F;p* < 20 GeV, the contribution of double-resolved events was increased and the 
contribution of direct events has been reduced, whereas the contribution of single-resolved 
events was left unchanged. The correction factor C{r) in each F;^* bin is calculated according 
to the relative contribution of each process determined in this way. 


The F;p* dependence of the jet shapes in di-jet production is presented in Figure The 
predictions of PYTHIA for direct, single-resolved and double-resolved processes and their sum 
are compared to the measured jet shapes. The central value of the plotted data points is 
obtained with the correction calculated from PYTHIA. The systematic error was obtained from 
the variation of the EGAL energy by ±5% and from the difference of the results from PYTHIA 
and PHOJET. The relative contribution of the processes in PYTHIA has also been changed 
over a range consistent with what the match to Figure will allow, and the resulting changes 
in the corrections have been added to the systematic error. The jets become narrower as F^p* 
increases. The predictions of PYTHIA with the default relative contributions of the different 
processes using the SaS-lD parametrisation reproduce the data reasonably well. There is almost 
no difference between the predictions of PYTHIA and PHOJET using GRV (not shown). The 
differences of the jet shapes for direct, single- and double-resolved processes are expected to be 
due to the different fractions of quark and gluon jets. It has been shown that gluon jets are 
broader than quark jets According to the prediction of PYTHIA using SaS-lD the fraction 
of gluon jets is 12% for direct events, 16% for single-resolved events and 41% for double-resolved 
events. 


The fraction of the transverse energy of the jets inside a cone of radius r = 0.5 around 
the jet axis, = 0.5) is shown as a function of Fp* (Figure ^) and as a function of 
(Figure Qd). The position of the data points is the mean value of Fp* in each bin. In Figure ^ 
the data are compared to PYTHIA with and without multiple interactions and to PYTHIA 
without initial (ISR) and hnal state QGD radiation (FSR). The PHOJET prediction is also 
shown. The difference between PYTHIA with and without multiple interactions (mi) is very 
small, whereas the PYTHIA prediction without initial and hnal state QGD radiation, where 
only fragmentation effects have been taken into account, leads to jets which are signihcantly 
narrower. This difference increases with increasing Fp*. Fig. shows that there is no observed 
dependence of the jet shape on xhe prediction of PYTHIA is in good agreement with the 

data. Gluon jets are predicted to be broader than quark jets. The hadron jets in the Monte 
Garlo events have been identihed as a quark or gluon jet depending on the type (quark or 
gluon) of the closest parton in the rjcj) plane. 

The results for '0(r) are presented in Figure ^ and d for both x^ smaller and larger than 0.8. 
It is observed that the measured jet shapes for x^ < 0.8, where more gluon jets are expected, 
are broader than those for > 0.8. For both regions of x^ the Monte Garlo generators 
reproduce the data reasonably well. 
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9 Inclusive two-jet cross-sections and NLO calculations 


To obtain absolute jet cross-sections which can be compared to theoretical calculations, we use 
the Monte Carlo simulation and an unfolding program to correct for the selection cuts, 
the resolution effects of the detector and the background from non-signal processes. The same 
technique was used as described in Ref. |]l|. To improve the performance of the unfolding 
program in the region > 3 GeV, bin-to-bin migration effects from jets at low must be 
taken into account. Therefore the jets are actually found with = 2 GeV and the unfolding 
is performed in the full > 2 GeV range, but the unfolded jet cross-section are only shown 
for > 3 GeV. The uncorrected number of jets reconstructed in the range 3 < < 4 GeV 

is about 15% larger with = 3 GeV than with = 2 GeV. This difference decreases to 
less than 5% for higher T^p*. 


In Fig. ^ the inclusive two-jet cross-section is shown as a function of i^p*. The average 
transverse energy, (F'p*), within each bin, which is plotted on the abscissa, is determined as 
proposed in Ref. |^. It is obtained by integrating an exponential function which is fitted to 
the neighbouring data points. The error on (i?p*) is calculated by varying the slope of the 
exponential function. The error bars show the statistical and the systematic errors, calculated 
in the same way as in Section and an additional error from the unfolding procedure are added 
in quadrature. The bin sizes, which are indicated by the tic marks at the top of the Figures, 
approximately reflect the experimental resolution. The results are summarised in Table ^ 


The F^p* distribution is compared to a NLO perturbative QGD calculation of the inclusive 
two-jet cross-section by Kleinwort and Kramer [^] who use the NLO GRV parametrisation of 
the parton distribution functions of the photon Their calculation was repeated for the 
kinematic conditions of this analysis. The renormalisation and factorisation scales are chosen to 
be equal to i?p\ The scale dependence of the NLO QGD calculations is expected to be small 0. 
The strong coupling constant Og is calculated from the two-loop formula with = 130 MeV, 

iVio 

f5) 

since this value is also used in the NLO-GRV parametrisation. Ghanging from 130 to 250 

iVlo 

MeV only in the Og formula increases the two-jet cross-section by factors from 1.4 to 1.07 in 
the range 3 < E^r^^ < 16 GeV for = 130 — 136 GeV |2^ . 


The data points are in good agreement with the calculation except in the hrst bin where the 
calculation predicts a much higher cross-section. The symmetric cuts on i7p* lead to singulari¬ 
ties of the NLO calculations. This problem only affects the first bin, where the NLO calculations 
are not reliable. The NLO QGD calculation gives the jet cross-section for massless partons, 
whereas the experimental jet cross-sections are measured for hadrons. The uncertainties due 
to the modelling of the hadronisation process have not been taken into account. Because the 
partons in the Monte Garlo models and the partons in the NLO calculations are dehned in 
different ways it is impossible to use the Monte Garlo to correct the data so that it can be 
compared with the NLO parton level predictions. If PYTHIA had been used to calculate a 
correction we would have had to increase the cross-section by a factor of between 1.2 and 1.3, 
with the largest effects at low i7p*. 


The predictions for the direct, single- and double-resolved parts and their sum are shown 
separately. The resolved cross-sections is the largest component in the region ^ 8 GeV, 
whereas, at high E^^^ the direct cross-section is largest. 
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10 Influence of the underlying event 


The NLO QCD calculations also do not take into account the possibility of an underlying 
event which leads to an increased jet cross-section. The underlying event is simulated in the 
Monte Carlo models PYTHIA and PHOJET which will be used to compare to different LO 
parametrisations of the parton distribution, SaS-lD 0, GRV 01 and LACl |0]. These 


sets of parton distributions contain different parametrisations of the gluon density with LACl 
predicting a much larger gluon density than GRV and SaS-lD. In PYTHIA and PHOJET 
the modelling of the underlying event includes multiple interactions. The contribution from 
multiple interactions has to be tuned using quantities which are not directly correlated to the 
jets, since otherwise effects of the parton distributions and of the underlying event cannot be 
distinguished. A signihcant difference between the predicted two-jet cross-sections obtained 
with PYTHIA and PHOJET using the same parametrisation of the parton distributions was 
observed in our studies at lower energies |]^. By adjusting the cutoff parameter p™ for multiple 
interactions the model dependence should decrease signihcantly. 


It is expected that the transverse energy flow outside the jets measured as a function of 
is correlated to the underlying event 0]- No effect due to the underlying event is expected 
for direct events at large x^. The increase of the transverse energy flow outside the two jets at 
small x^ can therefore be used to tune the number of multiple interactions in the model. 


The events were boosted into their centre-of-mass system and the transverse energy flow 
was measured as a function of x^ in the central rapidity region |? 7 *| < 1. The regions around 
the jet axes with R < 1.3 are excluded from the energy sum. As in Section |^, x^ denotes 
that the transverse energy flow of each event is added to the plot at the values of xl^ and of 
x~. Fig. 0 shows the transverse energy flows corrected to the hadron level. The systematic 
error was obtained from the difference of the results from PYTHIA and PHOJET and from 
the variation of the ECAL energy by ±5%. Fig. 0a shows the results of PYTHIA using the 
LACl parametrisation with different cutoff parameters. The transverse energy flow for the 
default p™ of 1.4 GeV/c is much too high in the hrst x^ bins. Without multiple interactions the 
transverse energy flow is too low. An optimised value of p™ = 2.5 GeV/c leads to a reasonable 
description of the data. The average number of interactions in double-resolved di-jet events 
is 4.1 for LAGl with p™ = 1.4 GeV/c and 1.5 with p™ = 2.5 GeV/c. 


The difference of the transverse energy flow using SaS-lD in PYTHIA with and without 
multiple interactions (p™‘ = 1.4 GeV/c) is very small and the predicted transverse energy flow 
is in good agreement with the data. The best description using PYTHIA and GRV is obtained 
with a p™^ of 2.0 GeV/c. For all further comparisons with PYTHIA, the cutoff parameter p™‘ 
was set to 2.5 GeV/c for LAGl, to 2.0 GeV/c for GRV and to 1.4 GeV/c for SaS-lD. 


In 7 p collisions at HERA the GRV parametrisation with a cutoff parameter of p™ = 
1.2 GeV/c for PYTHIA has been found to be in good agreement with the data, whereas 
the cutoff parameter was set to 2.0 GeV/c for LAGl |^. With this choice the models slightly 
overestimate the transverse energy flows at low x^ in our data. 

PHOJET with either SaS-lD or GRV is in reasonable agreement with the data. Ghanging 
the default cutoff of PHOJET from p™‘ = 2.5 GeV/c does not affect the transverse energy flow 
signihcantly, the parameter p™ was therefore left unchanged in PHOJET. 
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11 Inclusive two-jet cross-sections as function of |? 7 j®^| 


The size and dependence of the inclusive two-jet cross-section, which is dominated by 
the low events, depend on the chosen parton distribution functions which mainly dif¬ 
fer in the assumptions on the gluonic content of the photon. This leads to different pre¬ 
dictions for the inclusive two-jet cross-section, especially for double-resolved events. The 
inclusive two-jet cross-section as a function of Ir/j®*! is shown in Figure H for events with 




> 4 GeV and > 3 GeV and in Figure || for events with > 5 GeV and 

> 3 GeV. The asymmetric cuts have been applied because NLO calculations with sym¬ 
metric cuts are not infrared safe 




31 


The data sample is separated into events with a large contribution from double-resolved 
processes by requiring < 0.8 (Fig. P? and |b) and into events with a large contribution from 
direct processes by requiring > 0.8 (Fig. ||c and|^c). The results are summarised in Tables |] 
and The average values are consistent with being at the centre of the bins. Each jet 

is included with its value of in the cross-section measurement. Within the statistical and 
systematic uncertainties of the measurement, the data distributions are nearly independent of 
pjg_ p with a small drop towards Ir/j®*! = 2, whereas the data distribution of the 

direct events shown in Fig. p; at = 2 falls to about half of its value at Ir/j®*! = 0. In the 
kinematic range shown, this is in agreement with the expectations of the Monte Garlo models. 


The NLO QGD calculation of the inclusive two-jet cross-section is in excellent agreement 
with the differential cross-section da/dlr/j®*! shown in Figs. |^a and ^a. The SaS-lD parametri- 
sation with the PYTHIA and PHOJET models predicts a two-jet cross-section which is 
signihcantly too low for the whole data sample and for events with a large contribution from 
double-resolved events 


^7 


xy < 0.8). The cross-section using GRV m is a bit too low for the 


whole data sample but is able to describe the cross-section for the events with x^ < 0.8. In this 


range the cross-sections as predicted using LAGl are much too high. It should be noted 


that the overall normalisation of jet cross-sections can alway be adjusted in a LO calculation 
by changing the value of the strong coupling Og. 


In contrast to these discrepancies between the predicted cross-sections for x^ < 0.8 due 
to different parametrisations of the parton distributions, the differences between the predicted 
cross-sections for events with a large contribution of direct events (x^ > 0.8) seem to depend 
mainly on the model, PHOJET or PYTHIA, and not the parametrisation used (Figs, p: and p:). 
The PHOJET cross-section is slightly too high, whereas PYTHIA describes the data well. As 
expected, this is independent of the chosen cutoff parameter for multiple interactions. 


To further reduce the influence from multiple interactions and hadronisation effects the 
inclusive two-jet cross-section was also measured for events with E^^^ > 5 GeV. (Fig. |l^). The 
results are summarised in Table The fraction of resolved events is smaller in this data sample 
and the difference between the different parametrisations decreases. Nevertheless the LAGl 
prediction is still too high for the whole data sample and especially for the data sample with a 
large contribution of double-resolved events (Fig. |T^b). The fraction of the cross-section with 
x^ > 0.8 compared to the total jet cross-section is signihcantly larger in the range > 5 GeV 
than for the lower E^^^ ranges. 
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12 Conclusions 


We have measured di-jet production in photon-photon interactions with the OPAL detector 
at e^e“ centre-of-mass energies of 161 and 172 GeV with an integrated luminosity of 
20 pb“^. Jets were identihed using a cone jet hnding algorithm with i? = 1 in the kinematic 
range > 3 GeV and < 2. Two-jet events originating mainly from direct and double- 

resolved photon interactions were separated experimentally using the variables and x~. The 
Monte Garlo models PYTHIA and PHOJET describe the transverse energy flow around the 
jets reasonably well. 


The distribution of the parton scattering angle 9* has been reconstructed from the rapidities 
of the two jets. Data samples with large direct and double-resolved contributions have been 
compared to LO and NLO QGD calculations. A strong rise has been observed in the | cos 0*1 
distribution of the data sample with a large double-resolved contribution at high | cos 0* |, as ex¬ 
pected from QGD. The flatter | cos0*| distribution of the data sample with a large contribution 
from direct events is also in good agreement with the QGD calculation. 

The energy prohle of the jets has been measured in different regions of . It has been 
observed that jets with high are narrower than jets with small values of E^r^^. The jet shape 
in events with x^ < 0.8, where double-resolved events are expected to dominate, is found to 
be broader than the jet shape in events with > 0.8, where direct events are expected to 
dominate. These differences are assumed to be caused by different fractions of quark jets and 
gluon jets. The influence of multiple interactions and of different parametrisations of the parton 
distribution functions of the photon is small. 


The inclusive two-jet cross-sections were measured as a function of and |? 7 p*|. The 
measured cross-sections are in good agreement with next-to-leading order QGD calculations by 
Klasen, Kleinwort and Kramer |T^ above E^^^ = 4 GeV using the NLO GRV parametrisation 
of the parton distributions of the photon. 


The inclusive two-jet cross-section is dominated by the resolved processes in the low 
region. In order to distinguish between the contributions to the jet cross-section from possible 
multiple interactions between the spectator partons and from the parton densities, the contri¬ 
bution of multiple interactions in the models has been tuned using the measured transverse 
energy flow outside the jets at low x^. However, within the errors of the measurement we are 
unable to differentiate between models with and without multiple interactions. 


Using PYTHIA and PHOJET the LO GRV parametrisation is also able to describe the 
two-jet cross-section whereas the cross-section predicted based on the SaS-lD parametrisation 
is too low and the prediction based on the LAGl parametrisation is signihcantly too high. 
The same behaviour is observed using a data sample with a large contribution from double- 
resolved events. As expected, there is no sensitivity to the choice of parametrisation for the 
complementary data sample with a large contribution from direct events. This behaviour still 
holds if the inclusive two-jet cross-sections are measured for > 5 GeV. 
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cos 6** 

da/d cos6** [pb] 

resolved {x^ <0.8) 

direct {xfj > 0.8) 

0.000 - 0.106 

0.106 - 0.213 

0.213 - 0.319 

0.319 - 0.425 

0.425 - 0.531 

0.531 - 0.638 

0.638 - 0.744 

0.744 - 0.850 

6.2 ±1.7 ±0.3 

2.7 ±1.1 ±0.2 

9.7 ±2.1 ± 1.1 

10.3 ±2.2 ±2.3 

10.5 ±2.2 ±1.7 

22.7 ±3.3 ±2.4 

30.1 ±3.8 ±2.5 

51.9 ±5.0 ±3.6 

4.2 ± 1.4 ±0.2 

5.0 ± 1.5 ±0.4 

5.6 ± 1.6 ±0.4 

7.8 ± 1.9 ±0.1 

6.8 ± 1.8 ±0.5 

7.2 ± 1.8 ±0.4 

6.5± 1.7±0.1 

15.8 ±2.7 ±0.7 


Table 1: Differential two-jet cross-section as a function of \ cos6'*|. The cross section is shown 
for events with < 0.8 and for events with > 0.8. The first error is statistical and the 
second error is systematic. 


(GeV) 

{Elf^) (GeV) 

da/dE^^ [pb/GeV] 

3.0 - 4.0 

3.47 ±0.01 

163 ± 3 ±16 

4.0 - 5.0 

4.47 ±0.01 

73.6 ± 2.0 ± 8.7 

5.0 - 6.5 

5.69 ±0.01 

27.9 ± 0.9 ± 3.9 

6.5 - 8.5 

7.42 ±0.01 

11.5 ± 0.5 ± 1.4 

8.5 - 11.0 

9.64 ±0.01 

3.83 ± 0.26 ± 0.60 

11.0 - 15.0 

12.73 ±0.02 

1.30 ± 0.14 ± 0.35 

15.0 - 20.0 

17.16 ±0.02 

0.12 ± 0.04 ± 0.13 


Table 2: The inclusive two-jet cross-section as a function of . The first error is statistical 
and the second error is systematic. 
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da/d V®* [pb] 



no xj^ cut 

resolved <0.8) 

direct ( 


> 

0.8) 

0.0 - 

- 0.4 

118.7 

± 

3.7 

± 

12.8 

77.2 

±3.3 

± 

11.8 

15.3 

± 

1.2 

± 

2.5 

0.4- 

- 0.8 

114.4 

± 

3.6 

± 

10.7 

73.3 

±3.1 

± 

10.4 

15.1 

± 

1.2 

± 

2.4 

0.8 ~ 

- 1.2 

93.4 

± 

2.9 

± 

11.7 

57.2 

±2.5 

± 

10.9 

13.6 

± 

1.1 

± 

2.4 

1.2 - 

- 1.6 

95.1 

± 

3.5 

± 

11.7 

56.0 

±2.8 

± 

11.5 

10.9 

± 

1.3 

± 

2.2 

1.6 - 

- 2.0 

81.2 

± 

3.8 

± 

9.7 

60.5 

±3.4 

± 

10.9 

8.2 

± 

1.2 

± 

1.9 


Table 3: The inclusive two-jet cross-section as a function of\'rf^^\ for events with > 4 GeV 
and >3 GeV. The inclusive two-jet cross-section is shown for all two-jet events and for 
data samples with a large contribution from double-resolved events by requiring < 0.8 and 
with a large contribution from direct events by requiring > 0.8. The first error is statistical 
and the second error is systematic. 



da/d ?7j®* [pb] 

no xfj cut 

resolved <0.8) 

direct {xfj >0.8) 

0.0 - 0.4 

0.4 - 0.8 

0.8 - 1.2 

1.2 - 1.6 

1.6 - 2.0 

58.3 ±2.6 ±7.7 

57.8 ±2.4 ±7.6 

49.8 ±2.0 ±6.5 

52.8 ±2.4 ±6.6 

45.2 ±2.9 ±7.2 

35.0 ±2.3 ±6.5 

33.6 ±2.0 ±6.0 

31.3 ±2.0 ±6.7 

28.8 ±2.0 ±5.8 

30.2 ±2.3 ±6.2 

8.8 ±0.8 ±1.9 

9.5 ±0.9 ±1.8 

9.2 ±0.9 ± 1.9 

7.9 ±1.0 ±1.7 

5.9 ±1.0 ±1.5 


Table j: The inclusive two-jet cross-section as a function for events with > 5 GeV 

and >3 GeV. The inclusive two-jet cross-section is shown for all two-jet events and for 
data samples with a large contribution from double-resolved events by requiring x^ < 0.8 and 
with a large contribution from direct events by requiring x^ > 0.8. The first error is statistical 
and the second error is systematic. 


Kl 

da/d pj®* [pb] 

no xf^ cut 

resolved < 0.8) 

direct >0.8) 

0.0 - 0.4 

0.4 - 0.8 

0.8 - 1.2 

1.2 - 1.6 

1.6 - 2.0 

31.0 ± 1.4 ±3.5 

32.0 ± 1.6 ±3.4 

27.6 ± 1.4 ±3.2 

27.5 ± 1.6 ±5.8 

18.3 ± 1.5 ±3.1 

17.2 ± 1.1 ±4.4 

15.6 ± 1.1 ±2.9 

11.3 ±0.9 ±3.3 

11.0 ±0.9 ±3.1 

7.5 ±0.9 ±3.3 

8.3 ±0.9 ±1.7 

9.3 ±1.0 ±1.8 

9.2 ± 1.0 ± 1.7 

7.8 ±1.0 ±1.6 

5.8 ± 1.2 ± 1.5 


Table 5: The inclusive two-jet cross-section as a function of\ri^^^\ for events with Elf^ >5 GeV. 
The inclusive two-jet cross-section is shown for all two-jet events and for data samples with a 
large contribution from double-resolved events by requiring xfj < 0.8 and with a large contribu¬ 
tion from direct events by requiring xfj > 0.8. The first error is statistical and the second error 
is systematic. 
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Figure 1: Uncorrected distribution in bins of the mean value of where is calculated 
as the mean value of of the two jets with the highest Elf^. The background has been 
subtracted using the Monte Carlo. The data points are compared to the predictions of PYTHIA 
(continuous line) and PHOJET (dashed line). The hatched histogram is the direct contribution 
to the PYTHIA events. The Monte Carlo histograms are normalised to the number of data 
events. Statistical errors only are shown. 
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Figure 2: Uncorrected energy flow transverse to the beam direction measured relative to the 
direction of each jet in two-jet events and normalized to the total number of jets, Nj^t, in the 
sample. Jets from data samples with a large contribution of (a) double-resolved and (b) direct 
events according to their xfj and x~ values are shown. The energy flow is integrated over 
|A0| < n/2. Statistical errors only are shown. The data points are compared to the PHOJET 
(continuous line) and PYTHIA (dashed line) simulations. 
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|cos6*| 

Figure 3: Angular distribution of events with large direct and large double-resolved contributions 
according to the separation with and xf. The data are compared (a) to LO QCD matrix 
element calculations and (b) to NLO QCD calculations using the GRVparametrisation. The 
curves are normalised to the data in the first three bins. The open circles show the distribution 
of events with x^ > 0.8 and the full circles show the distribution of events with x^ < 0.8. The 
inner error bar shows the statistical error and the outer error bar the statistical and systematic 
errors added in guadrature. 
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r r 

Figure 4-' The measured jet shapes, 'ip{r), correeted to the hadron level for eaeh of the two highest 
jets. The jet shapes are shown in bins of (a) 3 < E^rf^ < 6 GeV, (b) 6 < E^^ < 9 
GeV, (e) 9 < E^^^ < 12 GeV and (d) 12 < E^rf^ < 20 GeV. The predietions of the direet, single- 
resolved and double-resolved proeesses and their sum as predieted by PYTHIA are shown. The 
inner error bar shows the statistieal error and the outer error bar the statistieal and systematie 
errors added in quadrature. 
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Figure 5: The fraction of the transverse energy of the jets inside a eone of radius r = 0.5 around 
the jet axis is shown (a) as a function of and (h) as a function 0 / 77 ^®*. The measured jet 
shapes corrected to the hadron level, ipir), are shown in (c)forxf^ < 0.8 and in (d)forxf^ > 0 . 8 . 
The statistical error is smaller than the symbol size. The error bars show the statistical and 
systematic errors added in guadrature. 
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Figure 6: The inclusive two-jet cross-section as a function of for events with < 2 
compared to the NLO calculation by Kleinwort and Kramer mi- The direct, single-resolved 
and double-resolved cross-sections and the sum (continuous line) are shown separately. The 
inner error bar shows the statistical error and the outer error bar the statistical and systematic 
errors added in guadrature. 
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Figure 1: (a) Transverse energy flow outside the jets in the central rapidity region |? 7 *| < 1 as 
a function of x^. The statistical error is smaller than the symbol size. The error bars show the 
statistical and systematic errors added in guadrature. 






















OPAL 



Figure 8: The inelusive two-jet eross-section as a funetion of\r]^^^\ for events with > A GeV 
and > 3 GeV are shown for (a) all events and (b) for events with a large eontribution 
of double-resolved events by reguiring < 0.8 and (e) for events with a large contribution of 
direct events by reguiring > 0 . 8 . The curve in (a) shows the prediction of the NLO QGD 
calculation using the NLO GRVparametrisation. The inner error bar shows the statistical error 
and the outer error bar the statistical and systematic errors added in guadrature. 
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Figure 9: The inelusive two-jet cross-section as a function of\r]^'^^\ for events with > h GeV 
and > 3 GeV are shown (a) for all events and (b) for events with a large contribution 
of double-resolved events by reguiring < 0.8 and (c) for events with a large contribution of 
direct events by reguiring > 0 . 8 . The curve in (a) shows the prediction of the NLO QGD 
calculation using the NLO GRVparametrisation. The inner error bar shows the statistical error 
and the outer error bar the statistical and systematic errors added in guadrature. 






































Figure 10: The inclusive two-jet cross-section as a function of\rj^‘^^\ for events with > b GeV 
are shown (a) for all events and (h) for events with a large contribution of double-resolved events 
by reguiring < 0.8 and (c) for events with a large contribution of direct events by reguiring 
xff > 0 . 8 . The inner error bar shows the statistical error and the outer error bar the statistical 
and systematic errors added in guadrature. 
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